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Abstract 

In order to cope with the rapidly increasing service demand in cellular networks, more cells are 
needed with better resource usage efficiency. This poses challenges for the network planning since 
service demand in practical networks is not geographically uniform and, to cope with the non-uniform 
service demand, network deployments are becoming increasingly irregular. This paper introduces a new 
idea to deal with the non-uniform network topology. Rather than capturing the network character (e.g. 
load distribution) by means of stochastic methods, the proposed novel approach aims at transforming 
the analysis from the physical (irregular) domain to a canonical/dual (uniform) domain that simplifies 
the work due to its symmetry. To carry out this task, physical and canonical domains are connected 
using the conformal (Schwarz-Christoffel) mapping, that makes the rich and mature theory of Complex 
Analysis available. The main contribution of this paper is to introduce and validate the usability of 
conformal mapping in the load coupling analysis of cellular networks. 

Index Terms 

Cellular networks, conformal mapping, network planning, canonical domains, load coupling. 

I. Introduction 

Due to evolution of mobile cellular networks different generations of multiple access technolo¬ 
gies (from 1G to 4G) have been introduced [iT|. While it is not yet clear whether new air interface 
will be standardized for 5G [|2j], it is evident that densities of cellular networks are constantly 
growing to cope with the rapidly increasing service demand. Yet, geographical distribution of the 
service demand is not uniform, and accordingly, networks are becoming increasingly irregular. 
Thus, besides technical challenges on the air interface efficiency, network planning, modeling, 
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and performance analysis are becoming increasingly difficult, yet important, research fields. 

From the perspective of network analysis, each cellular deployment is certainly ‘unique’ in 
the sense that network layout, propagation conditions, and service demand are not identical in 
two different cellular deployments. Thus, reliable planning and analysis of the network usually 
takes a lot of time and computational effort since analysis is strongly dependent on the network 
scale and number of users (service demand volume) under consideration. Consequently, detailed 
system level simulations are not very suitable for network design/planning, but are more useful 
in performance evaluations and Radio Resource Management (RRM) testing. On the other hand, 
analytical approaches such as stochastic geometry E j4j, are much less costly to carry out and 
they have gained popularity as an alternative way to study the cellular networks. However, in 
most cases restrictive assumptions and heavy approximations are required, and hence, accuracy- 
tractability represents the fundamental trade-off between simulation-based and analytical methods. 

Our contribution is in the field of cellular network modeling and analysis. We present a 
novel solution framework that can be used for planning, analysis, and optimization purposes. 
Instead of using stochastic methods to capture the irregular nature of the network, we transfer the 
analysis to a dual domain that simplifies the analysis due to its symmetry. Thus, the idea aims at 
introducing a simple dual domain that is connected to original domain by means of a spatial 
transformation, more precisely, a conformal mapping [Q. Let us consider the general problem of 
finding network deployment (base station locations and corresponding cell areas) when service 
demand distribution and total traffic volume are known in a certain spatial area, namely, physical 
domain. In the proposed approach we form a conformal (Schwarz-Christoffel) transformation 
from a given (polygonal) physical domain to a rectangular area, referred to as canonical domain. 
Finding the (deterministic) uniform network deployment in (periodic or non-periodic) canonical 
domain to cope with the given service volume is straightforward and can be carried out using 
methods known from literature. The deployment compatible with the service demand in the 
physical domain is then obtained by using the inverse of the original conformal mapping. 

While idea behind the proposed approach is simple, it introduces some heavy challenges: 
computation of conformal mappings and their corresponding inverse mappings fulfilling the 
requirements of our problem is not straightforward. Even more importantly, the selection of the 
mapping parameters such that, e.g., service demand distribution is preserved, is challenging and 
related to the concept of conformal equivalence. We emphasize that the proposed methodology, 
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although being very promising, will require in future more in-depth treatment by means of 
classical complex analysis 0. Finally, we note that to the best of the knowledge of the authors, 
the idea is novel and provides a new perspective for cellular deployment studies, i.e., a new line 
of research. Therefore we have put more emphasis on explaining the ideas rather than technically 
presenting all the mathematical details, particularly those related to numerical methods. 

The rest of the paper is organized as follows: the next section describes in more details the 
research problem and the proposed solution. Spatial transformations through conformal mapping 


are presented in Section III and in Section IV the canonical domain and the analysis therein is 
explained. Numerical examples and the corresponding explanations are provided in Section |Vj 


Finally, Section VI closes the paper with conclusions and discussion on the future work. 


II. Problem Statement and Proposed Solution 
A. General description 

Cellular networks provide wireless access for users that are usually nonuniformly distributed 
within the network service area, by using a spatial partition of the area into cells. Traditionally, 
the planning has been done so that the operator is able to provide a certain Quality of Service 
(QoS) in terms of data rates and (high) coverage probability, i.e., users should be able to access 
their services and applications (almost) anytime and anywhere. Thus, the planning process starts 
from the characterization of the service demand. Given its importance, service demand patterns 
have been extensively studied in the context of cellular networks 07] (8j. 

In [|9]], a notion of irregularity in terms of the compatibility between service demand and 
capacity provision was presented. In the ideal case, the network should absorb the service demand 
with maximum resource usage efficiency. If all serving base stations have the same amount of 
resources, then cell coverage areas should be planned such that cells are equally loaded. Since 
the demand is nonuniformly distributed, such ideal topology would also be irregular with cells 
of different sizes. Motivated by these simple but fundamental observations, the authors have 
addressed the following research problem: 

Research problem: Assume a statistical description of the service demand in the area A, 
e.g., in terms of traffic distribution 5 and volume V. Determine the spatial partition 
A = A\ U • • • U Al of the network service area into cells such that the service volume 
Vi in each cell Ai is the same. Then select the base station locations accordingly. 
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Fig. 1. Research framework: Obtaining a deterministic reference model of the network that 1) is specific for the service demand 
spatial distribution (S) and volume ( V ), and 2) can be used for planning, load analysis, and layout/topology optimization. 


While this problem can be addressed by various computational methods, we are proposing a 
novel approach based on the classical complex analysis. To that end, two domains are considered: 

• The physical domain , characterized by the area A with service demand distribution 5 and 
volume V, and 

• the canonical domain , characterized by the rectangular area TZ with uniform service demand 
distribution 5' and the same service volume V as in the physical domain. 

In addition we need a spatial transformation F : A —>■ TZ that preserves the load coupling 
analysis ifTOll in different domains. Natural choice is to accomplish the connection between the two 
domains by means of conformal mappings 0, or more precisely, by means of Schwarz-Christoffel 
transformations mum- The big picture is shown in Figure [I] Details of the different phases of 
the solution will be provided in the next sections. The idea can be summarized as follows: 

1) Assume a certain spatial service demand distribution 5 and volume V. Service demand 
patterns are well-known, or at least good estimations are available [J7] [8]]. 

2) Deduce a conformal mapping F : A ^ TZ such that resulting service demand distribution 
S' in 7 Z is uniform and service demand volume is preserved. 
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Fig. 2. Representation of the required conformal mapping F and its inverse F 1 


3) Create a network deployment in the canonical domain. That is, define the number of base 
stations according to the service demand volume V. 

4) Deduce the inverse F~ 1 : TZ —> A that is used to map the uniform base station grid from 
the canonical domain back onto the physical domain. 

The phases 1) and 3) above are straightforward, however, the most challenging task in this 
approach is to find a suitable conformal mapping F in phase 2) and its inverse F _1 in phase 4). 

Finally, we note that the resulting network topology (base station locations and cell areas) 
provides useful information that can be used for planning, indicating the locations of the access 
points, or as reference in topology/energy optimization, e.g., in cases where a hyper dense network 
is deployed in A. 


B. Problem formulation and notations 

As it was indicated, the problem at hand, that of determining the spatial mapping F is difficult 
because several conditions need to be simultaneously fulfilled. Figure [2] illustrates the required 
mapping indicating the complex variables z and w used in the physical and canonical domain, 
respectively. In general, the transformation F and its inverse F _1 are expressed as follows: 

F(z) = w = u(x,y) + iv(x,y) = £ + irj, (1) 

F~ 1 (w) = z = w _1 (£, rj) + i rj) — x + iy, (2) 

where F and F~ l are conformal mapping^] if u, v, u~ l , and v -1 are harmonic (analytic) 
functions fl6[. Thus, the problem can be formulated as follows: 

'The notion of conformality is explained in the appendix. 
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TABLE I 
Notations. 


Symbol 

Description 

Symbol 

Description 

A 

Canonical domain 

1Z 

Physical domain 

S 

Spatial service demand distribution 

S' 

Uniform service demand distribution 

V 

Service demand volume 

V, 

Service volume of the I th cell 


Session time 

A 

Inter-arrival time 

L 

Number of base stations 

P min 

Minimum rate per user 

Ai 

Area of the I th cell in A. 

TZ;, 

Area of the Z th cell in TZ. 

F 

Conformal mapping: A—rTZ. 

U, V 

Harmonic components of F 

F - 1 

Conformal mapping: 1Z — > A 

u , v 

Harmonic components of F~ 1 

Q 

Quadrilateral 

m(Q) 

Conformal module of Q 


9 

Mapping from rectangle onto strip: 

f 

Mapping from rectangle onto polygon: C 

)> 

H 

Height of the canonical domain 1Z 

W 

Widht of the canonical domain TZ 


A h 

Vertical shift in 7Z 

A w 

Horizontal shift in TZ 

dP 

Periodic distance in TZ 

R 

Radius of hexagons in TZ 

Pi 

Power received from the 1 th cell 

r 

A given point in TZ 

i 

Average SINR 

V 

Average spectral efficiency 

T .bs 

'l 

Location of the 1 th base station 

13 

Propagation exponent 


find F : A—> TZ and F 1 : TZ —» A, 
subject to: 


5:A^ 


with 


5 = 1, 


IA 


S'(w) = 1/m, Vw e 7 Z, 
[ 5 ds = I S'ds'. 


(3a) 

(3b) 

(3c) 

(3d) 


We note that in ( |3c| ), m refers to the aspect ratio of the canonical domain and S, S' are area 
elements in the physical and canonical domain, respectively, as illustrated in Figure [2j In general, 
it is not difficult to find numerically some mappings F and F~ 1 between A and 7 Z lH3ll . However, 
Constraints (3b), ( f3c| ), and ([3d]) make the task hard because not only any conformal mapping 
F needs to be found, but the mapping should make the spatial service demand distribution 5 


uniform {5') in the canonical domain TZ. In (3d), the integrations correspond to standard area 
integrals of the densities 5 and S' over S and S', respectively. Hence, S' is the image of S after 
applying F and S is the image of S' after applying F -1 , as shown in Figure [2] The rationale 
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for formulating the mapping of A onto a rectangle 1Z will be formally justified in Section IV 


Finally, in order to facilitate the reading of the rest of the article, the summary of the notations 
is presented in Table [TJ 


III. Conformal Mapping 

In general, given a function f(x) and another function k(x. s), it is always possible to define 
the transformation: F(s) = J x k(x. s)f(x)dx as long as the product k(x. s)f(x) is integrable on 
the set X. In this work, the use of conformal mapping, a certain spatial transformation in the 
complex plane, is motivated by the following important features: 

/ The existence of the transformation and its inverse is guaranteed. An important and profound 
result from complex analysis, the Riemann’s mapping theorem (see [0 p. 221]) guarantees 
the existence of a conformal mapping between between non-empty, open, and simply 
connected proper subsets of C, such as A and 7 Z in Figure [2j 

/ Conformal transformations preserve notions such as locality and proximity between domains, 
and preserve a certain structure that is required in our context. 

The theory of conformal mapping is rich and mature. Comprehensive treatments can be 
found in [|5l [6] [13]]. In any case, a brief introduction to the notion of conformality and mapping 
composition is provided in the appendix. Hereafter, the focus is on the type of conformal 
mapping used in this work: Schwarz-Christoffel transformations IfTTII . To keep the presentation 


comprehensive, we explain the basics of Schwarz-Christoffel transformations in Section III-A In 


Section III-B we focus on mappings from polygons onto rectangles. This discussion is essential 
from the practical perspective since it describes the tools used to carry out phases 2) and 4) of 


the research problem in Section II-A 


A. Schwarz-Christoffel transformations: essentials and practicalities 

Consider the mapping shown in Figure [3] Suppose that the vertices of the polygon Y are 
Wi, w 2 , ■ ■ ■ ■ v' n with the corresponding interior angles 0]7r, o 2 x. • • • , a n n and assume that these 
points map onto points zi, z 2 , ■ ■ - , z n on the real axis of the z- plane. Then a transformation F 
that maps the interior of Y onto the halp-plane TL + is obtained by using the Schwarz-Christoffel 
formula ifTTl : 
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Fig. 3. Schwarz-Christoffel mapping from a polygon to the half-plane (H + ). 


/ z n ~ 1 

U(z-z^dz, (4) 

k =i 

where A and C are constants that adjust the size and position of the polygon (T) generated 
by the mapping and the notation f~ indicates complex integration 0], i.e., integration over 
the complex variable z. The basic idea behind the Schwarz-Christoffel transformation (and its 
variations) is that the conformal transformation F may have a derivative that can be expressed as 
§zF(z) = n F k , for certain functions F k [|5] fill . A huge variety of conformal maps admit this 
model, and indeed, almost all conformal transformations with known closed-form expressions 
belong to this family. The use of Schwarz-Christoffel transformations involves three basic steps: 

1) The problem of determining the pre-vertices, zi, z 2 , • • • , z n , which are not known a-priori. 
This problem is also referred to as the Schwarz-Christoffel ‘parameter problem’ and it 
is a very challenging numerical task of any Schwarz-Christoffel mapping procedure. The 
pre-vertices influence on the lengths of the sides in T. Very often, the constants A and C 
are also included in the parameter problem. 

2) The second step involves the evaluation of the integrals in the parameter problem. Given 
that there are only few examples of Schwarz-Christoffel mapping that can be studied by 
means of closed-form expressions (typically n < 3), these analytical examples are best 
serving as a teaching tool. Fortunately, there are well-developed numerical methods |[5]|. 
We note that practical problems usually require the composition of several numerical 
Schwarz-Christoffel mappings, as in the case of interest herein. 

3) Finally, computational effort is needed for inverting the mapping, i.e., finding the inverse 
transformation. Indeed, finding the inverse mapping (F~ l ) is more challenging since no 
formula exists in general. The vast majority of the techniques proposed so far, are based on 
two main approaches: Newton iteration and numerical solution of the initial value problem. 





A complete description of the appraches can be found in IflTI and references thereirj^J 


B. Conformal equivalence among quadrilaterals and the polygon-to-rectangle mapping 


Consider the simply-connected domains A and A 1 shown in Figure 4a We note that A is 
simply-connected if a closed curve in A is the boundary of some area contained in A. A closed 
curve is a curve where the beginning and end points are the same. A closed curve is called 
linear if it is a concatenation of a finite number of segments, i.e., straight lines. A curve is 
said to be simple if it does not intersect itself. A simple closed linear curve is called polygon. 
Therefore, the domain A is bounded by a polygon while the domain A! is not bounded by a 
polygon, but a closed curve. A quadrilateral 0 is defined as a system composed by a domain 


(such as A or A') and four different points in its bounding curve. Thus, Figure 4a shows two 
different quadrilaterals, Q — { A: zi, z 2 , z 3 , z 4 } and O' = { z\ , zf z 3 , z' 4 }. The ‘conformal 
module ’ 0 rn{Q) of a quadrilateral 0 is defined as the unique value m (see Figure |4a| for 
which a conformal mapping between A and 1Z exists. The conformal module m(Q) is itself 
another unknown of the problem of determining F : A —> 1Z. We note that the aspect ratio of 
1Z is completely determined by 0, and hence, the conformal mapping is only possible onto a 
rectangle 1Z with a certain unique aspect ratio rn(Q). Two quadrilaterals Q and O' are conformally 
equivalent [0, if m(Q) = rn{()') lfT51l . Later on, the important role of m ( 0) will be shown when 
comparing the correspondence between the load coupling analysis in the canonical and physical 
domains. 

Figure |4b] illustrates the method proposed in ||T6ll , based on Schwarz-Christoffel transformations. 
It can be used for computing approximations to m(Q) and the mapping F~ l : TZ —>• A, in 
cases where the domain A is (bounded by) a polygon. The method is designed to avoid the 
crowding phenomenon [HU §2.6]. associated with the use of previous conventional methods]^} by 
applying an infinite strip (rather than the unit disc or the upper-half plane) as the intermediate 
domain. Thus, given a polygonal domain A, with vertices in counterclockwise order at the points 


2 This explains why the potential of Schwarz-Christoffel mapping was fully exploited much after its discovery, when computers 
became available. As a result, today there are several efficient numerical methods, and also a number of public domain software 
packages, such as 03, to assist researchers with the numerical computation burden that is required. 

’The required numerical computation of the conformal map is difficult, or impossible, when the unit disc or the half plane are 
used as intermediate domains m- 
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Quadrilaterals 




(a) Quadrilaterals and equivalence (b) Transformation composition: rectangle «->• polygon 

Fig. 4. Conformal equivalence among quadrilaterals is a fundamental requirement for polygon-to-rectangle mapping. 


Cl, C2, • • • • Qn 1 and interior angles a in, o 2 n. • • • , a n n, the method proposed in lfl6ll is based on 
the expression 

F~ 1 = fog:K^A, (5) 


where g : 1Z S conformally maps 1Z onto the infinite strip S = {z \m(z) e (0,1)}, and 
/ : S —)• A is a modified Schwarz-Christoffel transformation. The function g maps a rectangle 
onto a strip of width 1, sending the corners 0, 1, 1 + im(Q), and i m(Q) of the rectangle to the 
points i, 0, R, and R + i, respectively. Then, according to lfT6ll . 


z = g(w) = —log [sn(z|m) ]. 


7T 


( 6 ) 


In ([6]), sn(;:|m) corresponds to the Jacobi’s elliptic function [fl3l . which is defined in terms of 
the incomplete elliptic integral of the first kind as follows: 


sn(/c|m) = sin(^), 

with 

(7) 

rv 

dt 


' I y/r 

— m sin 2 (t) 

(8) 


Clearly, 0 in Figure |4bj is linked algebraically to the conformal module, so algebraic spacing of 
the points on the rectangle corresponds to algebraic spacing of the points on the strip. In this 
case, the conformal module is g and it is an unknown in the problem. Furthermore, the images 
of rectangle corners are constrained to lie on the vertical lines: Re(z) = 0 and Re(z) = /?. At 
this point of the overall process, there are n — 3 unknowns, to be determined by the n — 3 side 
length conditions as in the standard half-plane or disc map iflTll . The solution of the parameter 
problem implicitly determines the correct value of R, and hence, the elliptic parameter m = e 2nR 
and the conformal module m(Q). An appropriate formulation for solving the parameter problem 
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is given in Ifl6ll . Accordingly, the function / in (|5]) is given by 

f* Z ^ 

C = f(z) = A + C J JJ [sinh (y - z k ) 


i (“fe-i) 


dy, 


(9) 


k =1 


where the z k s are the pre-images on Im(^) = 0 and \m(z) = 1, of the wf s, k — 1, 2 • • • , n. 
Note that ([9]) is very similar to ([3]), except that the hyperbolic sines are included. 

Thus, the composition indicated in ([5]) is used for finding the mapping F in the second step 
of the procedure given in Section |II-A[ This approach solves problem ([3])-([3d]) in cases where 
the spatial service demand distribution can be expressed or approximated in terms of analytic 
functions u 1 and v~ 1 that compose the mapping F~ x , see The mathematical relationship 
between F and 5 is given by 

1 1 


6(z) =8(F-\w )) = K 


V • F 


-l 


= K 


+ + v -i (£ + i v y 


V z e A, (10) 


dtj 


where V- is the classic divergence operator and K is a constant to guarantee that (3b) is fulfilled. 
This implies that 5 is a real analytic function. Given that < |TT)[ ) is the inverse of the divergence, it 
can be understood as a sort of convergence operator, which is intuitively aligned with the fact 
that service demand is assumed to be concentrated on the area where the grid is tighter, see 
Figure [2j Finally, it should be noted that the mapping F and its inverse F~ 1 convey information 
of the spatial service demand distribution 8. This guarantees that the allocation of the regular 
network obtained in the canonical domain (discussed in the next section) is mapped onto the 
physical domain (through F ~ 1 ) consistently with the service demand, and more importantly, in a 
deterministic manner. 


IV. Analysis in the Canonical Domain 

The analysis in the canonical domain corresponds to the third step of the procedure given in 
Section II-A In what follows, we apply a load-coupling model [TTOTl to estimate the resulting 
uniform load in the canonical domain. However, other models such as full load , related to the 
worst-case scenario, can be used as well. As indicated before, the target here is to create a regular 
network with a number of base stations able to satisfy the service demand volume V. We note 
that, this part of the analysis is independent of the spatial service demand distribution 8 in the 
physical domain. 
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A. Requirements and description 

In m, the notion of irregularity is presented as a measure for the dispersion of average loads 
at different cells. The rationale for this concept stems from the observation that in an ideal 
scenario, where the demand is uniformly distributed, a regular network topology where each cell 
takes exactly the same service demand and receives the same amount of intercell interference, 
would result in a uniformly loaded network. Hence, minimizing the irregularity, as indicated 
in m, would imply that for a nonuniform demand distribution, an irregular network topology is 
needed to balance the load at each cell. 

In order to guarantee the same load level at each cell in TZ, meeting all following conditions 
suffices: 

1) The service demand share is the same for each cell, i.e., the coverage area of each cell 
is equal and the demand is uniformly distributed. Hence, if there are L base stations, the 
coverage of each cell is equal to \7Zi\ = where \7Zi\ is the size of the Z th cell in 7 Z. 

2) The distribution / r d of the Euclidian distances in 7 Z is the same for all the points r e 1Z. 

The previous conditions guarantee that each cell will receive exactly the same amount of 

intercell interference. Intuitively, in order to fulfill the first requirement, the radio access network 
needs to have a regular geometry, hence, any regular tiling, e.g., rectangles or hexagons, should 
be used. To meet the second requirement, spatial periodicity suffices. The flat torus provides an 
example of 2D manifolcj^] that is topologically equivalent to the rectangular domain JZ iflTIl . 

As it was indicated, the analysis in the canonical domain (7 Z) aims at finding the number L of 
cells that is required to satisfy the uniformly distributed service demand with volume V. The 
analysis presented herein focuses on hexagonal tiling. The formulation for others regular tilings, 
such as rectangles, can be easily inferred from the formulation presented herein. 

A representation of the canonical domain (JZ) is shown in Figure [5] As cells are uniformly 
distributed in JZ, uniform density equal to is obtained, where W and H are the width and 
height of the rectangle, respectively. The radius of each polygon is R, and cells are obtained by 
means of the corresponding Voronoi tessellation |(T8l that uses bases stations as generators. Thus, 
given L access points, the I th cell (JZi C JZ) is defined as follows: 

A[ = {r G JZ : ||r-n|| < ||r-rj, Vj G £ = { 1, 2, • • • , L }, j ^ l}. (11) 

4 A manifold is a topological space that resembles Euclidian space near each point. 2D manifolds are also called surfaces. 
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2D periodicity in the torus 
R: Radius An example with 12 cells 

Fig. 5. The canonical domain corresponds to the rectangular torus in which a regular tiling can be embedded. The torus, 
topologically equivalent to the rectangle, illustrates the notion of spatial periodicity, i.e., wrap-around. 


In general, the allocation of an arbitrary number ( L ) of polygons in 7 Z is not possible. 
Furthermore, it is required to compute the relative angles for the selected lattice (rectangular 
or hexagonal), which is not a trivial task ifTTTl . In addition, as mentioned before, 7 Z needs to 
fulfill a certain aspect ratio (^), which cannot be controlled for some regular tilings, such as the 
hexagonal. Fortunately, these constraints do not represent a major issue from a practical point 
of view, since small deviations from the theoretical number of base stations may not imply a 
significant performance variation. 

Now R = and obviously, \7Z\ —W-H. Although the periodicity guarantees that the 

load analysis in TZ is invariant to translation, keeping the base stations as far from the boundary 
of TZ as possible is convenient for practical reasons. The shifts A w and A h , see Figure [5j are 


given by ^ and respectively. Finally, after some elementary calculations, the periodic 
(shortest) distance between two points rr and r 2 e 7J was found to be: 


(W 



2 7r" 

IV 

f H 



2 7T 1 

/(^ acos 

cos 

( r w _ r w) 

W . 

J) + 

— acos 

^2 7T 

COS 

(r? - r\) 

' ~w\ 


( 12 ) 


The letters ‘w’ and ‘h’ refer the horizontal and vertical coordinates of r, r w and r h , respectively. 
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B. Load analysis 

In load analysis, the target is to obtain a relationship between the network density (proportional to 
L ) and the resulting uniform/average load (a c ) in the canonical domain 7 Z. It is assumed that there 
is one single type of service. Each flow is satisfied if a minimum (average) throughput p min can be 
provided. In general, the demand is distributed according to some spatial distribution S : X —>■ M, 
where X is the spatial domain under consideration. Hence, 

I 5(x) dx = 1. (13) 

J x 

Obviously, in case of the 1Z, S' is a uniform distribution, see Figure [lj The service demand 
volume (V) is modeled in terms of two random variables: inter-arrival (A) and session time (//). 
In this manner, V is expressed in terms of the average number of users in the system as 
follows: V = Recall that V is the same for both A and TZ, see Figure jlj 

In general, the average Signal-to-Interference plus Noise Ratio (SINR) at r e TZi can be 
expressed as follows: 

R* (r) 

li{r) = - --, (14) 


K l&C\{l*} 


Oil Pi(r ) +o 2 


where Pi(r ) is the average received power from the Z th base station and a 2 is the noise power. The 
index /* indicates the serving base station. The instantaneous load ( 07 ) at any given cell is defined 
as the fraction of the available bandwidth that is being used. In case of Orthogonal Frequency 
Division Multiple Access (OFDMA), ai is proportional to the ICI that is being generated by 
the cell. As it can be inferred from ( p~4] >, the intercell interference coming from a neighbor cell 
is proportional to its average load (on = E{o:/}). In many practical cases, the SINR can be 
approximated by the Signal-to-Interference Ratio (SIR), i.e., 

Mr) 


7 i( r ) 


(15) 


a iM) 

The average spectral efficiency (fj) at some point r e TZi, is expressed by means of a 
nondecreasing function (funk) of the SINR, such that fj(r) = fuak{'Ji(r) ) : AT —> M [^] . Thus, 
the bandwidth requirement of a user located in r e TZ t is given by bjr) = /) niin ///(r) [Hz]. 

The average load (ai) in the I th cell is now given by 

1 


on = 


B, 


V) b h 


(16) 


sys 






where, 


Vi = 


5(x) dx 


and 


b, = 


'K, 


5(x) 
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\ EM 

J E{A}’ 

(17) 

b u (x) dx [Hz], 

(18) 


In, \J ni Sc(w)dw, 

Here V/ and bi correspond to the average number of users and bandwidth consumption in the 
I th base station, respectively, and £> sys is the system bandwidth. In order to estimate the average 
load vector (a = [ di a 2 • • • oil}), the iterative algorithm proposed in [[9J can be used. Yet, 
given the conditions of 7Z, it is clear that all the cells have the same load o c . Furthermore, the 
terms P/r) in ([Td]) and (|~i~5j) depend on the periodic distance (gP), in this case between the point 


of interest (r) and the location of the I th base station (r bs ), according to (12). Thus, under the 
assumption that all the access points transmit the same power, ( p3| ) can be written as follows: 

Mr) 


7 i( r ) 


(19) 


a c(P T (r) Pi*( r )) ac ^ [d p (r, r bs ) ] ^ — [d p (r, r bs ) ] ^ 

where P T (r) is the total received power at r, Pi*(r) is the received power from the serving 
base station, and (3 > 2 is the propagation exponent. The location of the serving base station 
(r bs ) depends on the selection of IZi, see Figure [5} but in any case, <7 p (r, rjT) 6 (0, /?|. Since 7 Z' 
satisfies the conditions given in Section |IV-A| the complete statistic of the SINR can be obtained 
by analyzing the points within any of the L fundamental parallelograms (one of them is indicated 
in dashed-red in Figure [5]) by means of the formulation just given. Alternatively, analyzing 
the coverage of one single cell also suffices. Thus, a c can be obtained by analyzing only the 
points r G TZi. Since the sum of the distances, and hence the received power, are independent of 
the any translation of the grid of access points, it is possible to write 



2 7T 

COS ( rf w 


H 

+ ( — acos 

Z 7T 


i 2?r 
cos | d h — 


( 20 ) 


where it is assumed that L = L\y x Lh- The terms d w and d h in ( |20| ) are indeed functions of the 
w, h, and r. They are defined as follows: 

3 


d w (w,r) = r w — I (w — 1) - R + A v 


( 21 ) 


= r 


a/3 


— ( (h — 1) a/3 R + mod(in + 1, 2) — R + A 


d h (w, h,r ) 


( 22 ) 
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Altogether, assuming funk(x) = log 2 (l + x), the relationship between the resulting load a c and 
the network density (or equivalently L ) will be given by the following expression: 

nin E{/i} 1 /' 1 


L = 


B sys E{A} \1Z\ ■ a c 


lo g 2 [l + K. 


p i* O) 


where 


Tli 02 V Oc-[.Fr(r)-Pi+(r)] 

refers to rounding up to the closets feasible integer. 


dr 


(23) 


V. Numerical Examples 

To illustrate the proposed methodology, some numerical examples are provided. First, we 
introduce two example domains, related service demand distributions, and discuss the computation 


of the required transformation. In Section V-B we carry out the required load analysis in the 


canonical domain and in Section V-C we discuss the correspondence between the canonical and 
physical domains in the light of examples. 


A. Conformal mapping 

In order to illustrate aspects of interest, two different physical domains have been considered 
for numerical evaluations. The domains A\ and *A 2 are shown in Figures [6a] and [6bj respectively. 
As it can be seen, both domains correspond to simple linear curves, i.e., polygons, although 
with different number of vertices. Moreover, nonuniform spatial service demand distributions 
(<$! and 8f) are defined in each case. The Cumulative Distribution Function (CDF) of demand 
probability (considering 500 x 500 area elements in each case) is also shown in Figure [6c] 

Numerical composite mappings [[TT| . as described in Section III-B[ are used to redistribute 
these spatial service demand onto canonical domains (rectangles). Mappings are obtained by using 
the Schwarz-Christoffel software provided in lfT4ll . Spatial representations of applied mappings are 
illustrated in Figures [7a] and [7b] The vertices and prevertices of the composite Schwarz-Christoffel 
transformation are indicated in the boundary of each domain by small pink circles. The 4 corners 
of each quadrilateral (and their images) are indicated by larger red and green circles and squares. 
Note that, as it was remarked earlier, the conformal module, m(Q), is an important unknown to be 
found. In other words, the aspect ratio of the rectangle representing the canonical domain cannot 
be selected arbitrarily once the quadrilateral is fixed. However, the selection of the vertices and 
quadrilateral’s corners in the physical domain gives the degrees of freedom to find an appropriate 












16 



-2-10 1 2 3 


(a) Physical domain 1: At (b) Physical domain 2: A 2 

Fig. 6. Physical domains considered for numerical examples. 



Demand probability ( p ) 

(c) CDFs: spatial service demand 



(a) Conformal mapping Ai 1Z 1 (b) Conformal mapping A 2 IZ 2 

Fig. 7. Representation of the Schwarz-Christoffel mappings used as examples. 


mapping. This task needs to be carried out numerically and no exact method exists. In case of a 
rectangular tiling, matching the aspect ratio of 1Z and m ( 0) can be done easily, while hexagonal 
lattices have their own limitations IfTTl and ratio matching is not always possible. 
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Fig. 8. Load in the canonical domain (hexagonal tiling). Setting: \TZ\ = 6.84 x 4.90 « 33.58, E{/i} = 120 s, E{A} = 0.05 s, 
L = [64, 100, 144, 196, 256, 324, 400], B sys = 5.0 MHz, r min = 100kbps. 


B. Load analysis in the canonical domain 

Figure [8] illustrates the load in the canonical domain (1Z) for different propagation exponents (d). 
The vertical axis indicates the resulting uniform load a c as a function of number of base stations L, 
shown in the horizontal axis. Note that the vertical axis is given in logarithmic scale. The figure 
shows that, increasing the propagation exponent results in less ICI, and hence, lower load levels, 
i.e., more capacity. Results of Figure [8] are obtained by solving (23) for the setting indicated in 
the caption text. In general, the resulting load a c depends on the service demand volume 0, the 
aspect ratio of 7Z (that is a constraint from the conformal mapping), and the type the of regular 
tiling that is used fT71 . However, ( [23] ) always allows for establishing the relation among the 
variables of interest (a c , V, and L ) for any rectangular domain fulfilling the conditions given 
in Section |IV-A[ and hence, it is a gateway to get a relationship similar to the one shown in 
Figure [8] for any arbitrary physical domain through its corresponding canonical counterpart. 


C. Correspondence between the canonical and physical domains 

The analysis in the canonical domain allows for calculating the number of base stations ( L ) 
required to fulfill the service demand volume with a uniform load a c . Indeed, such analysis 
depends exclusively on V and it could be done in 1Z without computing F. However, computing 
F allows obtaining F -1 , and hence, here is when conformal mapping plays a key role. As 
indicated before, the idea is use the inverse mapping of F (that makes S regular), i.e., F~ 1 , to 
place the base stations in the physical domain given that F captures the irregular nature of the 






















18 



f = 0.718 
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No linear 
transformation 



(a) Site mapping: Ai —t 7Zi: hexagonal tiling (L = 36) (b) Site mapping: Ai —> IZi: rectangular tiling (L = 36) 




(c) Site mapping: Ai —» 1Zi\ hexagonal tiling (L = 180) 

Fig. 9. Mapping site’s locations from the canonical domain to the physical domain using Schwarz-Christoffel transformations. 


spatial service demand distribution, see Figure [I] The site mappings 7Z) Ai and 1Z 2 O A 2 
of our two examples for L = 36 and L = 195, respectively, are illustrated in Figure [9] Thus, 
Figures [9a] and [9c] show the resulting topology if a hexagonal tiling is employed. As it can be 
seen, the aspect ratio ^ of TZ cannot be fit with the conformal modules m ( 0 ) when regular 
hexagons are used, and hence, an intermediate linear transformation is required. The resulting 


rectangular tilings are shown in Figures 9b and 9d The figures also indicate the conformal 
modules, m(Q), of the physical domains, and aspect ratios, of the canonical domains. 

The corresponding load coupling analysis is shown in Figure [TOj Therein, Figures UOaHlOdl 
indicate the resulting load level of each cell in the canonical domain, with and without periodicity, 
and in the physical domain, without periodicity. It should be noted that, in all the cases, 
a c > max {a;}. V/ G [1, L\ where 07 refers to load in cell l when the load coupling analysis is 
carried out in the canonical domain without periodicity (the dashed orange pattern). We note 
that, the load a c in the periodic canonical domain represents the worst case scenario in terms of 
intercell interference that flows over the domain boundary. In nonperiodic canonical and physical 
domains, the load in the cells close to domain boundary becomes smaller, which is also shown 
in the oscillating load pattern (dashed orange and green curves). The main interest is placed 
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f . TZ\, periodicity 

- 1Z\, no periodicity 

- Ai, periodicity 




"] 
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(a) A\-1Z\. hexagonal tiling 


(b) rectangular tiling 


(c) A 2 -T 12 - hexagonal tiling (dj A2-7Z-2'- rectangular tiling 



(e) Spatial pattern: *4i@hex (f) Spatial pattern: „4i@rec (g) Spatial pattern: _42@hex (h) Spatial pattern: ^ 2 @rec 

Fig. 10. Numerical examples of load-coupling analysis: comparison of results between canonical and physical domains. 


on the correlation between the load coupling analysis (without periodicity) in both domains: 
1Z and A. That is, the correlation between the dashed-orange and solid-green patterns. It has 
been observed that, the better the match between the aspect ratio oflZ and the conformal module 
of A, the higher the correlation between both load patterns. This has been verified in numerous 
numerical examples, which obviously cannot be presented herein due to space constraints. Thus, 
it can be concluded that, the existence of the aforementioned correlation between the load pattern 
in both domains makes evident the consistency of the proposed method. The load coupling 
correspondence between the domains is achieved because the spatial transformation preserves well 
the relative service demand share of each cell. In Figures [TOb] and lOd where the fit is perfect 
the correlation is almost perfect cell-by-cell. The periodic pattern is related with the fact that 


boundary cells receive, in general, less intercell interference as it is illustrated in Figures lOe lOh 


Furthermore, in cases where the aspect ratio of 1Z is clearly different from the conformal module 
of A, see Figures [9a] |10a| and [9c] |10c| the distances between average load levels in both domains 
are larger, although preserving the same periodic pattern cell-by-cell. Small differences are caused 
by the fact that the inverse conformal mapping ( F ~ l ) does not guarantees perfect service demand 
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sharing once the corresponding Voronoi regions are calculated, but good approximations of 
them. Another source of relative error, in case of certain regular polygons, such as hexagons, 
comes from the fact that no periodicity can be assumed in the physical domain, and accordingly 
the service demand share of the boundary cells varies more significantly as it can be seen in 


Figure lOg This problem is obviously alleviated when rectangles are used as it is shown in 


Figure lOh From a practical point of view, the fact that a c is (provided that conformal module 
is preserved in the canonical domain) above the actual load levels in the physical domain is quite 
advisable because it provides a natural protection against instantaneous variations with respect to 
average service demand values. In addition, if large scale networks are analyzed by parts, in a 
clustered fashion, it also provides a sort of wrap-around to account with the ICI coming from 
neighbor clusters. The resulting network topologies, which are deterministic in nature, can be 
used for planning and optimization of new or existing deployments, respectively. 

Finally, it is worth saying that, the resulting dispersion in the load pattern (central cells tend to 
be more loaded) due to the lack of interference coming from neighbouring cluster can be reduced 
by increasing the noise power arbitrarily in cells on the area boundary. With such adjustment, 
the correlation will not be altered and the periodic patterns (orange and green) would tend to 
match with the red line, i.e., a c . 


VI. Conclusions and Future Work 

This paper introduced a framework based on the idea of analyzing cellular networks through 
dual domains and spatial (conformal) transformations. The idea is novel and it provides a new 
approach to design, model, and analyze cellular network topologies in a deterministic manner, 
which is of great interest from both theoretical and practical point of view. Moreover, such 
network topologies are function of the service demand that has been characterized for the physical 
domain of interest, i.e., on the target coverage area. In this way, the service demand is modeled 
statistically, while the physical domain (the network coverage area) is represented deterministically 
by means of polygons. 

The numerical examples provide evidence that conformal mapping can be used as a tool to 
generate the required spatial transformations, while the load coupling analysis verifies a certain 
equivalence between the network topologies in the canonical and physical domain. Key conditions 
that must be fulfilled in order to improve the accuracy of the method have been also identified. 
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Given that the work in the canonical domain is, in general, easier than the analysis of irregular 
topologies, the framework presented herein opens the road towards a novel analysis. That is, 
the presented idea basically opens a new research framework, and consequently, the authors 
acknowledge that this is just the beginning and there is a plenty of future work ahead with 
several interesting open challenges. Among the most interesting ones are to postulate and verify 
some performance bounds, which are quite challenging in this context because some of the 
mathematical structures can only be dealt with through numerical methods. In addition, composite 
tessellations that can be used to model heterogeneous networks is, without a doubt, another 
interesting, yet challenging, research line that could be covered through the umbrella of idea 
presented herein. Finally, the rich framework of potential theory and boundary value problems 
will be of great help to understand or even device new spatial mappings that can also be used in 
the analysis of cellular networks through canonical domains. Given that conformal mapping and 
boundary value problems are strongly connected and have been used in other practical contexts, 
such as physics, this work also opens a promising line of research to formalize the notion of 
network irregularity, as the conformal mapping F conveys the information about the structure 
and density of the network topology that is created. Current research efforts are aligned in that 
direction. 
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Appendix 

Conformality and mapping composition 
A. The notion of conformality 

A function g : M n —* M n is called conformal if it preserves angles. What exactly does it mean? 
In the Euclidean norm, the angle between two vectors is defined by their dot product. However, 
most analytic maps are nonlinear, and so will not map vectors to vectors since they will typically 
map straight lines to curves. However, if we interpret ‘angle’ to mean the angle between two 
curves, or more precisely, the angle between their tangent vectors at the point of intersection as 
illustrated in Figure 1 la[ then it is possible to make sense of the conformality notion. Thus, in 
order to realize complex functions as conformal maps, it is needed to understand their effect on 
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(a) conformal transformation (b) Complex curve and tangent (c) The Riemann mapping theorem. 

Fig. 11. A conformal transformation establishes a unique (one-to-one) mapping from points ao £ A to other points a' Q € A', 
where the function / is analytic. The points where / is not analytic are called critical points. 


curves. In general, a curve in the complex plane is parametrized by a complex-valued function that 
depends on a real parameter t as follows: z(t) = x(t) + i y(t), t e [a, b]. In particular, the tangent 
vector to the curve at the point z can be identified with the complex number §f — §f + i fjr- 
Smoothness of the curve is guaranteed by the requirement that 7 / 0. When we interpret 
the curve as the motion of a particle in the complex plane, so that z{t) is the position of the 
particle at time t, the tangent ff 7 ^ 0 represents its instantaneous velocity. The modulus of the 


tangent, ||f| = y (jjj) + (||) , indicates the particle’s speed, while its phase, Z.z, measures 
the direction of motion, as prescribed by the angle that the curve makes with the horizontal, see 
Figure [TTb| The (signed) angle between two curves, C\ and C 2 , is defined as the angle between 
their tangents at the point of intersection z = z\{ty) — z 2 (t 2 ). If 0\ = Zzi(ti) and 0 2 = Zz 2 (t 2 ), 
then the angle 9 at z is 6 = d 2 — 6\. Now, consider the effect of an analytic map f(z). A curve 
C parametrized by z(t) will be mapped to a new curve D = f(C), where D(t) = f(z(t)). The 
tangent to the image curve is related to that of the original curve by the chain rule: c ^f i = 
Therefore, the effect of the analytic map on the tangent vector % is to multiply it by the complex 


dt 

number If the analytic map satisfies our key assumption 7 ^ 0 , then || 7 ^ 0 , and so the 


image curve is guaranteed to be smooth. Accordingly, 


I d/O) I _ I d/Q) dz |_| df(z) I I dz 


dt 


dz dt I 


the speed of motion along the new curve f(z(t)) is multiplied by a factor 


dz 1 1 

I d/O) 

I dz 


f |. Hence, 


> 0. It is 


worth noting that the magnification factor 


I d/011 

I dz I 


depends only upon the point z and not how 


the curve passes through it. Indeed, all curves passing through the point z are speeded up (or 


slowed down if 


I d/O) I 

I dz I 


< 1) by the same factor. Similarly, the angle that the new curve makes 

’^ 7 ^ + Z||. Therefore, the tangent angle 
i.e., its tangent has been rotated through 


'd/O) _ / d/O) dz _ 

- dt ^ dz dt 

— 7® 


with the horizontal is given by 
of the curve is increased by an amount Y , )z 
angle 4>. Again, the increase in tangent angle only depends on the point z (and obviously, the 



























23 


transformation /), hence, all curves passing through z are rotated by the same amount o. This 
immediately implies that the angle between any two curves is preserved. 


B. Existence and composition of conformal mappings 

One of the hallmarks of conformal mapping is that one can assemble a large repertoire of 
complicated mappings by simply composing elementary mappings. This relies on the fact that the 
composition of two complex analytic functions is also complex analytic. Hence, the composition 
of two conformal maps is also conformal. However, a fundamental question is: Is there, in fact, 
a conformal map f(z) from an arbitrary domain A to the unit disk? The theoretical answer is 
the celebrated Riemann’s mapping theorem. If A is any simply connected open subset, not equal 
to the entire complex plane, then there exists a one-to-one complex analytic map f(z), satisfying 
the conformality condition ‘-jr f 0 for all z E A, that maps onto the unit disk B = \z'\ < 1. The 
proof can be found in []6l p. 221], Thus, the Riemann’s mapping theorem supports the notion of 


mapping composition as shown in Figure lie 
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